Tailoring wave nonlinearity through spatial composites 
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We propose and demonstrate theoretically a method to achieve and design optical nonlinear 
responses through a light-mediated spatial hybridization of different standard nonlinearities. The 
mechanism is based on the fact that optical propagation through a spatial composite of different 
nonlinear media is governed by an effective nonlinear response if the spatial scale of the sequence 
is much smaller than the light diffraction length. We apply our general approach to the significant 
case of centrosymmetric photorefractive crystal biased by a periodically modulated external voltage 
to predict strictly bending-free miniaturized soliton propagation. 

PACS numbers: 



Wave propagation affected by nonlinearity can produce 
a variety of interesting and useful effects, one of these 
being the formation of solitons. In general, these emerge 
as the product of a specific wave-matter interaction, by 
which intrinsic wave features, such as diffraction or dis- 
persion, are countered by the nonlinear response, to pro- 
duce particle-like dynamics. In this respect, solitons are 
a manifestation of a specific physical context, which can 
generally be associated with the nonlinear model that 
characterizes the wave-matter interaction. 

We here focus on the possibility of reversing this per- 
spective, i.e., that of how to produce a prescribed soli- 
ton phenomenon, or, equivalently, of how to devise and 
tailor an appropriate nonlinearity. A possible solution 
would be to change the nature of the material, tailoring 
its physical and chemical properties, but this in general 
involves a lengthy material development and does not 
provide versatility. We here propose a method to shape, 
reconfigure, and design a nonlinear response, specialized 
for the generation of solitons, based on the macroscopic 
spatial composition of standard material nonlinearities. 
In particular, in conditions in which this sequencing has 
a scale that is smaller than the characteristic beam prop- 
agation scale (e.g. the diffraction length), the resulting 
optical propagation is driven by an effective nonlinearity 
that, although a product of the underlying combination, 
can have different quantitative and qualitative features, 
i.e., a natural hybridization produced by the system wave 
dynamics. 

The physical underpinnings are in many respects an 
elaboration of those that in mechanics are associated with 
the motion of a mass particle subject to the combined 
action of a time-independent force and a rapidly varying 
external one, a situation first investigated by Kapitza 
P, Here, even though the particle is not able to 
follow the rapid force oscillations, yet these produce an 



additional effective slowly-varying force field. The con- 
cept has been successfully extended and implemented in 
the general context of nonlinear physics [s*], such as for 
kink solitons in perturbed sine-Gordon models, with ap- 
plications for long Josephson junctions guiding-center 
solitons in perturbed Nonlinear Schroedinger Equations, 
with applications in optical pulse propagation in fiber , 
in quadratic nonlinearities with applications for spatial 
quadratic solitons ^ and in "nonlinearity management" 
in both Bose-Einstein Condensates dynamics and non- 
linear optics Q- In all these cases, the mechanism is 
based on a perturbation of the model propagation equa- 
tion either through an external modulated source term, 
or through a modulation of a linear and/or nonlinear 
medium property. What we propose here is to bring this 
Kapitza mechanism even further, by combining and se- 
quencing, at a fast scale, different media characterized 
by their own nonlinearities, thus achieving a qualitatively 
different nonlinear response (governing the average field) 
which is not a leading known nonlinearity and corrections 
due to modulations. 

Our proposal can find implementation in a vast variety 
of systems where macroscopic spatial sequencing is fea- 
sible. For example, a Kerr layer can be interlaced with a 
saturable system, such as a liquid- crystal cell, or simply 
with a Kerr layer with different nonlinear parameters. An 
example of a natural setting for this nonlinear designing 
is found in photorefractive crystals This is because 
the nonlinear response is not fixed by the nature of the 
material itself, but is also the result of the geometry of 
the applied external bias voltage. For example, in some 
conditions the optical nonlinearity can be switched from 
a self-focusing saturated response to a self-defocusing one 
simply by changing the orientation of the field [^l . In dis- 
tinction to standard nonlinear materials, this malleability 
allows spatial sequencing in a single crystal. 
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Consider the case of the scalar propagation of a parax- 
ial beam through a medium that has the nonlinear re- 
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sponse 

Sn{I,Ix,z) = cq{I,Ix) + ^ CmiI,Ix) COS f — — j , 

m— 1 ^ ^ / 

(1) 

where z is the propagation direction, / = is the 

optical intensity, = dl/dx and L„ the longitudi- 
nal scale of variation. Here we have modelled the lon- 
gitudinal sequence along the z-axis of different nonlin- 
ear responses by means of the superposition of a back- 
ground mean nonlinearity cq and a periodic nonlinear 
response whose period is Ly. Note that Eq.((T|) depends 
on the spatial derivative Ix in order to encompass also 
slightly nonlocal nonlinear responses. In the situation 
where the optical wave is not Bragg-matched with the 
periodic response Tc?], the backward reflected field can 
be neglected so that the slowly-varying part of the ID 
optical field A{x, z) satisfies the parabolic wave equa- 
tion [id/dz+{l/2k)d'^ldx'^]A = -{k/no)Sn{I,Ix,z)A 
where k — 27rno/A (wave- vector carrier), A the wave- 
length, no is the sample background index of refraction. 
Assume now that Ly is such that \ <^ Ly <^ Ld (the first 
inequality assuring the paraxial description) so that, re- 
tracing the spirit of the Kapitza mechanical approach 
P, [1], it is possible to set A{x, z) — Ao{x, z) + 5A{x, z), 
where Aq is that part of the field that has a longitudinal 
scale of variation L^, and 5 A is i) longitudinally rapidly 
varying, at a scale Ly, and ii) it is uniformly in the condi- 
tion \5A\ <C \Aq\. This self-consistently amounts to a de- 
composition of the field into a slowly varying mean-field 
component and a rapidly oscillating and small correction, 
a direct consequence of condition Ly <C Ld- Substitut- 
ing this field requirement into the parabolic equation and 
exploiting the smallness of \5A\ as well as the large differ- 
ence between the two scales Ly and Ld of variation, it is 
possible to derive an equation for 5 A which can be solved 
and substituted into the equation for A^. The resulting 
nonlinear evolution equation for the " average" field Aq is 
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Ao is thus governed by a slowly- varying (i.e., at scales 
^ Ly) effective nonlinear response Sneff which is the 
superposition of cq and a correction consisting in a non- 
trivial function of cq and Cm contributions. Therefore, 
by appropriately choosing the sequence of longitudinal 
known nonlinearities, we are in the position to gener- 
ate effective nonlinearities which are qualitatively dif- 
ferent from their underlying constituents. The validity 




FIG. 1: Hybridizing a nonlinearity: functional dependence of 
index change on the intensity, where Sua — Ual /{I + la) and 
Sni, = ni,I / {I + It) are two standard saturable nonlinearities, 
and various effective nonlinear responses Sueff — Co for differ- 
ent spatial compositions. The result is that different portions 
of a single wave experience simultaneously profoundly differ- 
ent self-action, i.e., focusing, defocusing, and even linear (i.e. 
none). Here we have chosen Ua — 5 ■ 10~*, la = 0.05(/a -|- h), 
nt = -1.3 • 10"^ and h = 0.95(/a + h)- 



of the present approach is guaranteed by the conditions 
Ly -C Ld and {Ly/X)J2m=i'^rn/m <C f (resulting from 
the requirement \6A\ ^ |^o|) the last inequality being 
accessible since, even though Ly A, the nonlinear coef- 
ficients Cm are generally small. Note that, in Eq.Q, the 
correction to cq is proportional to Ly, so, as we would 
expect, more rapidly oscillating patterns produces lesser 
corrective terms to the slowly-varying part of the nonlin- 
ear response. 

We are now in a position to fully appreciate the versa- 
tility of the idea. Consider the case of two different sat- 
urable nonlinear media, characterized by the responses 
6na = nal/{l -)- la) and Sni, = ml / {I + h) [na and n^ 
being the maximum saturated refractive index changes 
and la and lb the saturation intensities), sequenced along 
the z— axis with the period Ly. From Eq.([3]) we get 
Sneff — Co = [fSna + (f — f)Sn(,], where / is the frac- 
tion (0 < / < I) of the period occupied by material a. 
Profiles of Sn^ff are reported in Fig. (I) for a particular 
choice of the parameters of the two underlying saturable 
nonlinearities. Note the non trivial behavior of Sn^ff, 
being simultaneously self-focusing, self-defocusing, and 
even linear for different intensities, i.e., for different parts 
of the beam: a response that, to the best of our knowl- 
edge, does not correspond to any known material optical 
response, and whose features (maximum position, width, 
etc.) can be easily tuned by an appropriate choice of 
constituents saturable media (i.e. na, la, ni, and /;,) and 
of macroscopic /. 

In order to pit this method on a pressing issue of some 
import to recent experiments, we use it to devise an ac- 
cessible scheme to achieve photorefractive solitons that 
do not self-bend [HI, In fact, self-bending impedes 
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FIG. 2: Devising a bending-free self-trapping nonlinearity in photorefractives, (la), (2a) and (3a): Crystal layers, electrode 
geometries (gray and black stripes) and optical beam (yellow regions) configurations, (lb), (2b) and (3b) : Steady state optical 
intensity profiles. (Ic), (2c) and (3c) Nonlinear refractive index profiles supporting the corresponding optical propagations. 
Note the transition into the Kapitza-like regime in case (3). 



such important effects as miniaturized self-trapping p^ . 
i.e., the formation of solitons with micron-sized widths 
up to the ultimate paraxial Umit, and Umits counter- 
propagating soUtons, i.e., head-on collisions or vector 
soliton composites [IJ]. Moreover, bending can have 
maj or impact on the interaction of solitons in patterns 
[l5| , because it alters nonlinearly the relative orientation 
of the beam and the pattern. Consider a centrosym- 
metric photorefractive sample as in the three geome- 
tries reported in Figs.2-(la), 2-(2a) and 2-(3a) where 
the external bias potential is delivered on the facets 
X = Lx and x = —L^ by means two standard plate 
electrodes 2-(la), or a system of alternating electrodes 
along the propagation z-axis 2-(2a), 2-(3a). In the 
(l-l-l)D configuration at steady state, the light experi- 
ences the nonlinear refractive index change Sn — aE^ 
where a = —(l/2)nQg€^{er — 1)^, g is the significant 
quadratic electro-optic coefficient and Ex{x, z) is the x— 
component of the electric field due to the photorefrac- 
tively stored charges and to the presence of the exter- 
nal bias potential. The associated electrostatic potential 
^(x, z) (defined by the relation E = — V$), in condi- 
tion where charge saturation can be neglected [l^, is 
governed by the equation V • [QV$ — X^Q] = where 
X = KbT / q, (5 = 1 + 1 /Id, I{x, z) is the optical intensity 
and Id is the intensity of a reference background uniform 
illumination. Assuming the external bias potential to be 
of the form $(±La;,z) — ±Vbcos(K^z), in the condition 
where the transverse beam width a (along the x- axis) 
is much smaller than ~ l-n jn,^, the equation for $ 



admits the approximate solution 



cosh(K„a;) 
cosh(K^Lj;) 



1 
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dx' 



Qix',z) 



cos{ktjz) 
(4) 

which, for Ky = reproduces the well-known one- 
dimensional photorefractive potential in the standard 
condition of 2-(la) Q. Evaluating the x-component of 
the electric field (retaining only the main contributions 
and exploiting again the condition a << Ly) we obtain 
the nonlinear refractive index change 



5n = a 



X 



dl 



I + Iddx 



■ COS{KyZ) 



(5) 



where ip = VoId/{Lx cosh{KyLx)), so that light dynamics 
is here governed by a nonlinear response belonging to the 
general class of Eq. ^ . Therefore we can apply the above 
general method to predict that, according to Eqs.([2]) and 
([3]), the optical field Aq, experiences the effective nonlin- 
ear refractive index change 



SUeff = 



2(/ + /<i)' 



[i;' + 2x'{Ix? 



(6) 



where we have approximated Sn^ff = cq (neglecting 
higher order contributions). Compared to the standard 
nonlinearity (i.e., for k„ — 0), our case has acquired a 
reflection symmetry (x — > —x) for an even intensity pro- 
file I(x) — /(— x), a symmetry that is intrinsically un- 
available for the standard photorefractive response. The 
absence of symmetry is the cause of self-bending for soli- 
tons of several microns, but for smaller ones, it actually 
inhibits their stable formation [13|, [IGj] . It is reasonable to 
conclude that the we now have the unique possibility of 
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exciting micron-sized photorefractive solitons, and that 
these will propagate in a straight line. 

To probe the uniform-to-fast-modulation-regime tran- 
sition, we have numerically integrated the full time- 
dependent photorefractive nonlinear optical model Q. 
In our numerical approach, at each instant of time, we 
calculate the electric field distribution induced by the 
boundary applied voltage and the photoinduced charge 
solving the {x, z) electro-static Poisson equation, and the 
corresponding optical field distribution determined by 
the electro-optic response through the parabolic equa- 
tion [iSj . We have chosen a crystal bulk (layer) of potas- 
sium lithium tantalate niobate (KLTN) (e^ = 3 • 10'*, 
g — 0.13r7i'*C~^, no = 2.4) at room temperature of thick- 
ness — 2 X 55 /im and length Lz — 1000 /im through 
which an initial Gaussian beam launched at z = with 
an half width at half maximum wq = 3.5 /im (with a 
diffraction length of ~ 250 /im) and a peak intensity 25 /d 
propagates. In the three electrode geometries of Figs. 2- 
(la), 2-(2a) and 2-(3a), the gray and black stripes are 
electrodes whose voltages are at 15 and —15 Volts, re- 
spectively: in 2- (2a) and 2- (3a) they are alternated along 
z with periods 600 ixm and 200 /im. In Figs. 2-(lb), 2- 
(2b) and 2-(3b) we report the intensity profiles |Ap of the 
evaluated optical fields at temporal steady state whereas 
in Figs. 2-(lc), 2-(2c) and 2-(3c) we plot the underlying 
refractive index patterns supporting the corresponding 
optical propagations. Configuration (1) corresponds to 
the standard soliton response where a 3.5 micron-sized 
optical beam experiences a bending deflection ~ 25 /im 
after ~ 4 diffraction lengths. In Figure 2-(lc) we re- 
port the corresponding refractive index, with its curved 
profile. Note that the beam bends toward the side of 
negative external potential and, most importantly, that 
reversing the applied voltage the resultant beam is the 
mirror image of the former one since the parabolic equa- 
tion with 5n of Eq.(I5]) (with = 0) is left invariant by 
the transformation ^ — > —tp and x ~x. When the ex- 
ternal applied voltage is slowly modulated along the z— 
axis (configuration (2)), the optical beam shows a global 
wiggling intensity profile. In fact, since the modulation 
is slow {Ly > Ld), the Kapitza mechanism is ineffective: 
the situation must be regarded as a cascading of differ- 
ent media, each with an almost uniform externally bias 



voltage producing its own asymmetric distortion, charac- 
terized by a snake-like response profile (see Fig.2-(2c)). 
It is in configuration (3), where the external applied volt- 
age is made to oscillate along the z— axis with a period 
very close to the diffraction length, that a wholly dif- 
ferent effect arises: the self-trapped beam forms along a 
straight line. The absence of beam deflection is the signa- 
ture that the Kapitza mechanism, for the nonlinearities, 
is in action. To underline this, note how the effective re- 
sponse is not a trivial product of a material response, but 
is the result of a detailed underlying spatially oscillating 
z-dependent pattern, as reported in Fig. 2- (3c), yet light 
propagation betrays only a negligible intensity oscillation 
(we chose a limiting case with L„ ^ Ld). As a final conflr- 
mation of the Kapitza picture, i.e. that configuration (3) 
is truly governed by the effective nonlinearity, we have 
integrated both the parabolic equation for A with the 
5n given by Eq. ([5]) and the parabolic equation of Eq. Q 
with 5nf,f f given by Eq. ^ with a gaussian boundary (at 
z ~ Q) field profile corresponding to that used in the sim- 
ulation of Fig. 2- (3a) 2- (3b) and 2- (3c), obtaining, in both 
cases, good agreement with the results of Fig.2-(3b). 

In conclusion, we have devised a general method to 
achieve an artificial or effective nonlinearity, having light 
propagate through a rapidly-varying (but macroscopi- 
cally accessible) spatial composite of media, each gov- 
erned by its own nonlinearity. The fact that the ef- 
fective nonlinear response can be profoundly different 
from the underlying nonlinearities of the composite al- 
lows an efficient design of nonlinearity and, in particular, 
the possibility of a feasible observation of a prescribed 
soliton manifestation. To underline versatility, we have 
analyzed briefly the composition of two saturable media. 
To implement the idea to overcome a present limitation 
in soliton investigation, we have considered photorefrac- 
tive propagation in a modulated external voltage regime, 
identifying the conditions for observing micron-sized soli- 
tons not suffering the standard deflection arising from the 
self-bending effect. This particular result is of fundamen- 
tal relevance also in applied photonics, since it at once 
brings us one step closer to ultimate circuit integration 
(up to non-paraxial scales) [l3|, and provides the means 
to achieve head-on soliton collisions, the avenue to self- 
splicing circuits. 
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